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We solve an one-dimensional stochastic model of interacting particles on a chain. Particles can 
have branching and coagulation reactions, they can also appear on an empty site and disappear 
Q „ spontaneously. 

This model which can be viewed as an epidemic model and/or as a generalization of the voter 
model, is treated analytically beyond the conventional solvable situations. With help of a suitably 

I* chosen string function, which is simply related to the density and the non-instantaneous two-point 

*vj ' correlation functions of the particles, exact expressions of the density and of the non-instantaneous 

two-point correlation functions, as well as the relaxation spectrum are obtained on a finite and 
i i. periodic lattice. 

O ■ PACS number(s): 02.50.-r, 02.50.Ey, 05.50.+q, 82.40.-g 

(U ■ 

Because of their important role in the description of classical interacting many-particle non-equilibrium systems, 
+1> \ reaction-diffusion (RD) models have been extensively investigated in the last decade 0B|. In lower dimensions, 
they provide relevant examples of strongly correlated systems which cannot be correctly described by mean-field-like 
c/3 \ approaches. In this sense satisfying comprehension of RD models in lower dimensions would require exact solutions, 
which are scarce, even in one spatial dimension. In some cases, however, certain RD models are known to be solvable. 
These cases can essentially be classified into four categories: (i) models for which the equations of motion of correlation 
functions are closed [3]; (ii) the free-fermion models B (or systems which can be mapped onto the latter, see pLpf); 
some other (one-dimensional) RD models can be solved by the (iii) Matrix Ansatz method H and some other by (iv) 
the interparticle distribution function (IPDF) method |7H9]; first introduced for the study of the diffusion-coagulation 
model (and its variants). It has also to be mentioned that the solution of various one-dimensional RD models have 
been obtained from the diffusion-coagulation models via similarity transformations |2pC||. It has been established 
that the latter solvable situations correspond to free-fermion systems 12] . 

The purpose of this work is to present a generalization of the IPDF method and to apply this technique to solve an 
one-dimensional stochastic model which is not solvable using conventional methods. The model under consideration 
exhibits a massive spectrum, implying an exponential approach towards the steady-state. The expressions of the 
r-o density and non-instantaneous correlation functions are determined. 

t**»« ■ Consider a periodic lattice of L sites on which (classical) particles interact. Each site is either empty (denoted 

by the symbol 0) or occupied by a particle at most, say, of species A (hard-core interaction). When a particle and 
a vacancy are adjacent to each other, a branching reaction can take place and the particle A can give birth to an 
offspring (A0 — ► A A and 0A — ► AA) with rate T\q = Tq\; another possible reaction is the death of the particle 
(A0 — ► 00 and 0A — ► 00) with rate T™ = r§5 ■ When two particles are adjacent, they can coagulate (AA —» A0 
and AA — ► 0A) with rate T\1 = T^. In addition, when two vacancies are adjacent, a particle can appear (birth 
process, 00 — > A0 and 00 — > 0A) with rate Tqq = TqJ. The system described above can be viewed as an epidemic 
model where particles can spontaneously appear/disappear, have an offspring and coagulate. It can also be viewed 
as a generalization of the voter model B , where the presence/absence of particle is associated to an opinion (yes/no) 
and each site is associated to an human being. According to the dynamics of the model, each individual changes his 
opinion at a rate proportional to the opinion of his neighbours. 

A particle (vacancy) at each of the L— lattice sites corresponding to spin down (up), the master equation of the 
model can be rewritten as an imaginary-time Schrodinger equation for a quantum spin-chain problem: jk\P(t)} = 
—H\P(t)), where \P(t)) = Yl{ n \ P({ n }>t)\{n}} describes the state of the system at time t (the sum runs over the 
2 L configurations) and H is the stochastic Hamiltonian (non hermitian) expressed in a spin-^ representation as 
H = J2f=i H J-j+ii with 

- H jJ+ i = r?° {(1 - n J+1 )(a+ - nj ) + (1 - n,)(a+ +1 - n J+1 )} 
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+ ^lo {(1 - n 3+l ){aj + n.j - 1) + (1 - nj)(aj +1 + n j+1 - 1)} 

+ r ii { n j(^j+i - n 3+i) + n j+ i(crf - rij)} 

+ r io {%-t-iO/ + n j ~ 1) + nj{vj+i + %+i ~ !)} ) C 1 ) 

where the a^- are the usual Pauli matrices and n 3 = ^(1 — <r|) is the density operator at site j. We also define the 
"left vacuum" (x\ = 2{ n >({ n }l- The probability conservation yields (x\H = 0. 
For the model under consideration, the equation of evolution of the density, from an initial state |P(0)), is therefore: 

j t (n 3 ) = -ixMH^u + H jd+1 )e- m \P(0)) = 2A + B((n j+1 ) + (n^)) - 2C(n J -) + D ((n 3 n J+1 ) + (n^-i)) , (2) 

where A = r$, B = T\\ - r o °, C = T™ + T™ and D = r°° + r o ° - (TJJ + rg). When D = and B ^ C, for 
a translationally- invariant system with initial density of particles (rij(0)) = /?(0), the solution of (0) simply reads 

(rij(t)) — c B + (p(0) — C A B ) e -2 ( C-S ) f . However, when D ^ 0, it is clear from (|2|) that the equation of motion 

of the correlation functions of the model give rise to an open hierarchy P,p|, which is not in general solvable. In 
addition, the Hamiltonian (||) is not quadratic and cannot, in general (excepted when T\q = T™ and T\^ = TJq, 
see P] for a complete classification of free-fermion systems) , be casted into a free-fermion form. Furthermore, this 
model cannot be solved by the traditional IPDF method which is not applicable J7|-[| in the presence of the processes 
A0 — * 00 ; 0A — v 00 and in the absence of processes A0 — ► 0A ; 0A — ► A0 (the latter should occur with the same 
rate as the coagulation rates [Q-|j]). 

To our knowledge, no exact results have been obtained for the model under investigation beyond the D = and 
free-fermion cases. In order to obtain the exact expression of the density beyond the latter conventional cases , we 
generalize the IPDF method introducing the following string function (x + L > y > x > 1): 

S x , y {t) = ((a - bn x ){a - bn x+ i) . . . (a - bn y -i))(t), (3) 

where a and b are non- vanishing numbers. When a = b = 1, S x ,y(t) is the empty interval functions used in the 
traditional IPDF method [|7]-||. The idea to solve the model under consideration here (with certain restrictions for 
the reaction-rates) is to choose suitable a and b in order to have a closed equation of evolution for S x , y (t). This is 
achieved by imposing the following ratio between a and b 
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1 + -FTO > 1 ( 4 ) 



and for the following reaction rates: 
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According to (0) and with the rates (J5J) , for the model under consideration on a periodic lattice of L sites, we have 
(1 < x < y < x + L): 

TtSx-vtt) = I (&+!,»(*) + S x , y -i(t)) + f {S x -i,y(t) + S x . y+1 {t)) - jS x , y (t) -(y- x)SS x>v (t); (l<x<y<x + L) 
i t S x , x+L {t) = -LSS XtX+L (t) (6) 

& x , x (t) = 1, 

where a = 2(aC - bA),(3 = -2D/6,7 = 2(5 + C) - S and 5 = ^A > 0. The prescription S x<x (t) = 1 is obtained 
requiring that S x ^ x+ i(t) = a — b(n x (t)) and using (0). 

The subcase T}J = T$ implies a = j3 = B = D = and we recover (C* ^ 0) (n x (t)) = Q ~ Sa: fc x+l(t) = % + 

(K(0))-£) e - 2C *. 

Hereafter we focus on the more general situation where (0) are fulfilled with T}J ^ TJq, and thus a^0,/J^0 . 

It is useful to consider the auxiliary function lZ Xiy (t) = /j, x ~ v S Xt y(t), where we introduce the complex numbers 
H = -isgn(a)\ji\ 1/2 and q = i\a/3\ 1/2 ^ 0. Notice that, because of (§), < ^ < \. With help of (|), we obtain the 
equation of motion of lZ x . y (t): 

f t n XtV (t) = | J2 e =±i {Kx+e, y {t) + Tl x ,y +e {t)} ~ 7^x,j/(t) - (y - ^^.^(i); (l<a;<y<a; + L) 
£llx,x+L(t) - -LSTZ X , X+L {t) (7) 

^x,x(i) = 1, 



The stationary solution of p) is obtained with the Ansatz: TZ x , y (oo) — AL.J y _ x+UJ {2q/5) + Bi,Yy_ x+ul {2q/S), where 
J v (z) and Y v (z) are the usual Bessel functions of first and second kind, respectively, Al and Bl are constants to be 
determined. Inserting the expression of TS.^ 2,(00) into (pi), we obtain u> — "f/S. Taking into account the boundary 
conditions TZ x . x (t) = 72. X)X (oo) = 1 and 1Z x ^ x+ l(oo) = 0, we get 

7 = Y L+ , /s (2q/S) 

L J L+l/s {2q/5)Y l/s {2ql8)-Y L+l/s {2q/8)J l/s {2q/5) {) 

£ = J L+l/s (2q/S) 

L JL + ^/6(2q/S)Y rY/s (2q/S)-Y L+7/s (2q/S)J^ /s (2q/Sy { > 

which provides the stationary expression for the string function: 

Woo) = [J>- x (A L J y _ x+l/s {2q/5) + B L Y y - x+1 , s {2q/5)) (10) 

According to the definition of the string function, the density of particles at site x is given by (n x (t)) = a ~ x - fc x+l1 - > 
and therefore the explicit stationary density of particles reads 

Moo)) - a ~ 53 ^ +lM = i [a - n (X L J 1+l/s (2q/5) + B L Y 1+j/s (2q/5))] . (11) 



In order to solve the dynamical part of (0), we seek a solution of the form lZ x ^ y (t) — TZ x ^ y (oo) = X^a r y x e Xqt ■ 

A I ~,A i 

y— l,x ~ ' y+l,x 
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Thus the equation (J7|), for l<x<y<x + L gives rise to the following difference equation: r y _ 1 x + r y 

i+ r y,x+i = 2 ( 1+{v ~ x)S - A] r* >x . With the notation E = 3^-, this equation admits r$ >x = AJ y - x - E (2q/6) 
BY v _ x ^ E (2q/ $) as solution, where A, B and the spectrum {E} are determined from the boundary and the initial 
conditions. Indeed, the boundary conditions 7Z XtX (t) = 1 and 4z7Z,x,x+l = —L51Z xx+ l require respectively, r xx = 
AJ- E (2q/S) + BY_ E (2q/S) = and £a (q\ - SL) er^ 1 [AJ L - E (2q/S) + BY L _ E {2q'/S)] = 0, i.e., 

AJ- E {2q/6) +BY_ E (2q/S) =0 

AJ L - E (2q/S) + BY L _ E (2q/6) = 0, l±ZJ 

The only non-trivial solution of this system (for which A ^ and B ^ 0) requires 

J L . E (2q/S)Y_ E (2q/5) - J_ E (2q/5)Y L _ E (2q/8) = 0, (13) 

or equivalently in terms of Lommel function ]12| : 

R L -i,i- E (2i\q\/6) = 0, (14) 

Thus the relaxation spectrum of the string-function of the model is obtained as the zeroes of the Lommel function 
(O). The latter admits [L — 1) zeroes which are symmetrically distributed around ^ (which is also an eigenvalue 
if L is even) and have a degeneracy L. To obtain the complete set of L[L — 1) + 1 eigenvalues, i.e., the relaxation 
spectrum {Ei}, i = 1, . . . , L of the string-function (and not the spectrum of the Hamiltonian (p])), we have also to take 
into account the eigenvalue qX = L5, which follows directly from the boundary condition 4i1Z x . x +l = —LSTL x ,x-\-l- 

To our knowledge there are no explicit results on the zeroes of the Lommel function of imaginary arguments. In 
order to have more explicit information on the spectrum, we use the formal analogy, first noticed by Peschel et 
al. []9|, which exists between the problem under consideration and the energy spectrum of an electron in a finite 
one-dimensional crystal in an electric potential of strength £n (here £ = 1) |ll| . 

To compute explicitly {Ei}, i = 1, . . . , L, we take advantage of the following eigenvalue-problem: 

(E - n)F n = V (F n _! + F n+1 ) ;{l<n<L) 

F = F L = 0, (15j 

where F n = V E [J E - n (2V)J- E (2V) - (-l) n J E (2V)J„- E {2V)} are eigenfunctions. The eigenvalues of (|l|) are ob- 
tained as the zeroes of the following Lommel function: Rl-ia-e^V) = [O. Choosing V = ^f-, the problem of 
determining the relaxation spectrum is reformulated as that of solving the eigenvalue-problem (Uq) . The latter can be 
recasted into the following form: A4\!F)) = E\T)), where A4 is a (L — 1) x (L — 1) symmetric (in fact antihermitian) 



tridiagonal matrix and \J-)) is a (L 
of the matrix M. is the following: 



l)-components column- vector: \T)) = {F n —\ F-i 
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For small systems the (L 
{Ei} 



3,3± 



5+4V 2 ±V9+24V 2 +4f 7T 
2 



1) distinct eigenvalues {Ei} of (Eq) can be computed analytically. For L = 6, we have 
where we still have to take into account the additional eigenvalue qX = L5. 



For larger matrices we had to proceed numerically. Our analysis (based on the spectrum of large matrices, with 
L < 1000), shows that the spectrum {Ei} (and therefore {qX}) is real and symmetric around -| which is an eigenvalue 
when L is even. The other eigenvalues are not generally integers, but for the central part of the spectrum (when 
eigenvalues which are close of -|), the eigenvalues approach integer values. This is not the case at the extremities 
of the spectrum. In particular, the smallest eigenvalue E* — min Ei {Ei} is not an integer and depends on the size 
of the system: E* = cl > 1. However, for L > 1, 6i ^ 6^, and E* is a constant: E* — e^ > 1. For L = 6, 
we have the exact result e L=6 = 3 - ^ 5+4^+^24^+4^ with x < £l=( , < 3 _ i^g + ^Tjj. This expression can 
be considered as an excellent approximation to systems of size L» 1 and in particular for e^. As an illustration, 
for the case Tqq = 3/10, TJq = 1/2, TJ5 = 1 an d T™ = 1/3, with the expression above, we obtain (analytically) 
cl=6 = 1.0823337683. For larger systems (L — 10,25,40,1000), we obtain numerically (with an accuracy of 10 -10 ): 
eio = ^25 = £40 = £1000 = 1.0823337697. 

Therefore, the long-time dynamics (of large systems, with L ^> 1) is governed by the eigenvalue E* = e^ « 

3- 
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(17) 



The equation ( p7| ) provides the inverse of the relaxation-time of the system pa]. 

With the knowledge of the spectrum {Ei}, i = 1, . . . , L, the expression of the density at site x reads: 



KW) ~ KM) - ^^e- (B ' 5+7)t lY^ Ez (2q/5)J L - El (2q/S) - J 1 _ Ei (2q/S)Y L _ Ei (2q/S)] , 



(18) 



The coefficients A Ei are obtained from the initial condition (in the translationally-invariant situation, where S XtV {t) 
Sy- X (t)) according to: 



M = Yl [ AA_1 ] i : J {Jn-E j {2q/8)Y L _ Ej {2q/8) - Y n _ Ed (2q/S)J L _ Ej (2q/d)Y (5 n (0) - S n (oo)) ff 



(19) 



where J\f is an hermitian L x L matrix which entries read (for details, see (l3|): Afij — 

Et=i (Jn-E^q/S)Y L - El (2q/S) - Y n _ Ez {2q/S)J L - Ei {2q/S))* (J n _ E .{2q/S)Y L _ Ej {2q/5) - Y n _ Ej {2q/S)J L - Ej (2q/5)). 
The long-time behaviour of (n x (t)) is obtained in retaining in (|l8|) only the term Ei — E* . 



The result (18) can be extended to the computation of (n x (t)n Xo (0)). To do this it suffices to take n xo |P(0)) (instead 



of |P(0))) as initial state in fllq ) and thus to replace the coefficients A Ei ( J19| ) by those computed in considering 

({]li=x(« - 6nj(0))}»»o(0)) instead of Sy- x (0). 

In this work we have proposed a natural generalization of the IPDF method, to solve (with some restrictions on 
the reaction-rates) an one-dimensional reaction-diffusion model which can be viewed as an epidemic model and/or a 
generalization of the voter model and that could not be solved by previous approaches. 

On a finite and periodic lattice, we have obtained the exact expression of the steady-state, of the dynamical part 
of the density and of the non-instantaneous two-point correlation functions of the model under consideration, which 



exhibits a massive and real relaxation spectrum. This means that the steady-state (|11|) of the system is reached 
exponentially with a relaxation-time which is determined explicitly in (J17|). 

We thank Laurent Klinger for his computational assistance. The support of the Swiss National Fonds is gratefully 
acknowledged. 
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